Statistical properties of fracture in a random spring model

Using large scale numerical simulations we analyze the statistical properties of fracture in the two dimensional random spring model and compare it with its scalar counterpart: the random fuse model. We first consider the process of crack localization measuring the evolution of damage as the external load is raised. We find that, as in the fuse model, damage is initially uniform and localizes at peak load. Scaling laws for the damage density, fracture strength and avalanche distributions follow with slight variations the behavior observed in the random fuse model. We thus conclude that scalar models provide a faithful representation of the fracture properties of disordered systems.Comment: 12 pages, 17 figures, 1 gif figur

A Fast and Efficient Algorithm for Slater Determinant Updates in Quantum Monte Carlo Simulations

We present an efficient low-rank updating algorithm for updating the trial wavefunctions used in Quantum Monte Carlo (QMC) simulations. The algorithm is based on low-rank updating of the Slater determinants. In particular, the computational complexity of the algorithm is O(kN) during the k-th step compared with traditional algorithms that require O(N^2) computations, where N is the system size. For single determinant trial wavefunctions the new algorithm is faster than the traditional O(N^2) Sherman-Morrison algorithm for up to O(N) updates. For multideterminant configuration-interaction type trial wavefunctions of M+1 determinants, the new algorithm is significantly more efficient, saving both O(MN^2) work and O(MN^2) storage. The algorithm enables more accurate and significantly more efficient QMC calculations using configuration interaction type wavefunctions

Role of the sample thickness in planar crack propagation

We study the effect of the sample thickness in planar crack front propagation in a disordered elastic medium using the random fuse model. We employ different loading conditions and we test their stability with respect to crack growth. We show that the thickness induces characteristic lengths in the stress enhancement factor in front of the crack and in the stress transfer function parallel to the crack. This is reflected by a thickness-dependent crossover scale in the crack front morphology that goes from from multiscaling to self-affine with exponents, in agreement with line depinning models and experiments. Finally, we compute the distribution of crack avalanches, which is shown to depend on the thickness and the loading mode.Peer reviewe

Effect of Disorder and Notches on Crack Roughness

We analyze the effect of disorder and notches on crack roughness in two dimensions. Our simulation results based on large system sizes and extensive statistical sampling indicate that the crack surface exhibits a universal local roughness of $\zeta_{loc} = 0.71$ and is independent of the initial notch size and disorder in breaking thresholds. The global roughness exponent scales as $\zeta = 0.87$ and is also independent of material disorder. Furthermore, we note that the statistical distribution of crack profile height fluctuations is also independent of material disorder and is described by a Gaussian distribution, albeit deviations are observed in the tails.Comment: 6 pages, 6 figure

Crack roughness and avalanche precursors in the random fuse model

We analyze the scaling of the crack roughness and of avalanche precursors in the two dimensional random fuse model by numerical simulations, employing large system sizes and extensive sample averaging. We find that the crack roughness exhibits anomalous scaling, as recently observed in experiments. The roughness exponents ($\zeta$, $\zeta_{loc}$) and the global width distributions are found to be universal with respect to the lattice geometry. Failure is preceded by avalanche precursors whose distribution follows a power law up to a cutoff size. While the characteristic avalanche size scales as $s_0 \sim L^D$, with a universal fractal dimension $D$, the distribution exponent $\tau$ differs slightly for triangular and diamond lattices and, in both cases, it is larger than the mean-field (fiber bundle) value $\tau=5/2$

Fracture roughness in three-dimensional beam lattice systems

We study the scaling of three-dimensional crack roughness using large-scale beam lattice systems. Our results for prenotched samples indicate that the crack surface is statistically isotropic, with the implication that experimental findings of anisotropy of fracture surface roughness in directions parallel and perpendicular to crack propagation is not due to the scalar or vectorial elasticity of the model. In contrast to scalar fuse lattices, beam lattice systems do not exhibit anomalous scaling or an extra dependence of roughness on system size. The local and global roughness exponents (ζloc and ζ, respectively) are equal to each other, and the three-dimensional crack roughness exponent is estimated to be ζloc=ζ=0.48±0.03. This closely matches the roughness exponent observed outside the fracture process zone. The probability density distribution p[Δh(ℓ)] of the height differences Δh(ℓ)=[h(x+ℓ)−h(x)] of the crack profile follows a Gaussian distribution, in agreement with experimental results.Peer reviewe

Role of disorder in the size-scaling of material strength

We study the sample size dependence of the strength of disordered materials with a flaw, by numerical simulations of lattice models for fracture. We find a crossover between a regime controlled by the fluctuations due to disorder and another controlled by stress-concentrations, ruled by continuum fracture mechanics. The results are formulated in terms of a scaling law involving a statistical fracture process zone. Its existence and scaling properties are only revealed by sampling over many configurations of the disorder. The scaling law is in good agreement with experimental results obtained from notched paper samples.Comment: 4 pages 5 figure

Fracture Strength of Disordered Media: Universality, Interactions, and Tail Asymptotics

We study the asymptotic properties of fracture strength distributions of disordered elastic media by a combination of renormalization group, extreme value theory, and numerical simulation. We investigate the validity of the “weakest-link hypothesis” in the presence of realistic long-ranged interactions in the random fuse model. Numerical simulations indicate that the fracture strength is well-described by the Duxbury-Leath-Beale (DLB) distribution which is shown to flow asymptotically to the Gumbel distribution. We explore the relation between the extreme value distributions and the DLB-type asymptotic distributions and show that the universal extreme value forms may not be appropriate to describe the nonuniversal low-strength tail.Peer reviewe

An Efficient Block Circulant Preconditioner For Simulating Fracture Using Large Fuse Networks

{\it Critical slowing down} associated with the iterative solvers close to the critical point often hinders large-scale numerical simulation of fracture using discrete lattice networks. This paper presents a block circlant preconditioner for iterative solvers for the simulation of progressive fracture in disordered, quasi-brittle materials using large discrete lattice networks. The average computational cost of the present alorithm per iteration is $O(rs log s) + delops$, where the stiffness matrix ${\bf A}$ is partioned into $r$-by-$r$ blocks such that each block is an $s$-by-$s$ matrix, and $delops$ represents the operational count associated with solving a block-diagonal matrix with $r$-by-$r$ dense matrix blocks. This algorithm using the block circulant preconditioner is faster than the Fourier accelerated preconditioned conjugate gradient (PCG) algorithm, and alleviates the {\it critical slowing down} that is especially severe close to the critical point. Numerical results using random resistor networks substantiate the efficiency of the present algorithm.Comment: 16 pages including 2 figure

Morphology of two dimensional fracture surface

We consider the morphology of two dimensional cracks observed in experimental results obtained from paper samples and compare these results with the numerical simulations of the random fuse model (RFM). We demonstrate that the data obey multiscaling at small scales but cross over to self-affine scaling at larger scales. Next, we show that the roughness exponent of the random fuse model is recovered by a simpler model that produces a connected crack, while a directed crack yields a different result, close to a random walk. We discuss the multiscaling behavior of all these models.Comment: slightly revise